Ion induced quark-gluon implosion 
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We investigate nuclear fragmentation in the central proton-nucleus and nucleus - nucleus collisions 
at the energies of CERN LHC. Within the semiclassical approximation we argue that each nucleon 
is stripped in the average process off "soft" partons and fragments into a collection of leading quarks 
and gluons with large p t because of the fast increase with energy of the cross sections of soft and 
hard interactions. Valence quarks and gluons are streaming in the opposite directions when viewed 
in the cm. of the produced system. The resulting pattern of the fragmentation of the colliding 
nuclei leads to an implosion of the quark and gluon constituents of the nuclei. The non-equilibrium 
state produced at the initial stage in the nucleus fragmentation region is estimated to have densities 
> 50 GeV/fm 3 at the LHC energies and probably > 10 GeV/fm 3 at BNL RHIC. 

■ One of challenging theoretical phenomena is the fact that perturbative QCD interactions become strong at suf- 
ficiently small x due to fast increase with energy of gluon densities. A number of models were suggested for the 

■ theoretical description of this new QCD regime see, e.g. [1-3] for reviews. A common feature of all these models is 
that in the case of proton scattering off heavy nuclei with large radius Ra, the high-energy cross section for most 
important quark-gluon configurations in the projectile nucleon reaches the black body limit(BBL) er tot = 2irR A [4]. 

I '— ^ ' They have this feature because, for the majority of the quark-gluon components of the nucleon wave function, the 
interaction cross section satisfies the condition: a2pR^ 3> 1. Therefore, most components are absorbed by the nuclear 
target. Here p ra 0.16 fm~ 3 is the nuclear density. For the central nucleus-nucleus collisions, the cross section of hard 
| interaction is even more enhanced by the factor [gA(x,Q 2 )/irR A ] /(gN/^rf^) ra A 1 / 3 (gA/AgN) ra 6 for A ra 200. Part 
of the enhancement disappears because of the leading twist nuclear shadowing phenomenon. Evidently, the soft QCD 

■^j- \ interactions in the nucleon-heavy nucleus collisions at the CERN Large Hadron Collider (LHC) should also be close 

OS ■ to the BBL because of cr tot (pN) ra lOOmfr. 

The numerical estimates indicate that already for x ~ 10~ 4 the perturbative cross section of the interaction of qq 
dipoles of transverse sizes d >0.15fm with nuclei calculated within the framework of the QCD factorization theorem 
reaches BBL while for the octet dipoles of similar sizes this regime may start at x about 10 times larger. For the recent 
numerical studies see [8]. Hence, the proton (ion) - heavy ion collisions at LHC will be qualitatively different from 

■ those at fixed target energies and as well as at BNL Relativistic Heavy Ion Collider (RHIC) where typical x > 10~ 2 
'■pj \ for the central rapidities. This should lead to a number of novel soft and hard phenomena [9,10]. 

Here we will argue that in this new regime, at the initial stage of the heavy ion collisions, a very dense non- 
equilibrium quark-gluon state will be formed in the nucleus fragmentation region. Before addressing the case of heavy 
ion collisions, we briefly discuss the case of the central pA collisions. To visualize the space time picture of the BBL 
in the proton-nucleus collisions, let us consider the rest frame of the nucleus. Within the semiclassical approximation 
, valence partons keep the momentum of energetic parent projectiles. This property of the infinite momentum frame 
k>( ' Schrodinger wave function is well known in the non-relativistic quantum mechanics, in the pQCD approximations, 
rS | in the parton model approach. It is valid also in BBL where a strong PQCD interaction is important for the large 
intervals in rapidity only. But because of dominance in this case of the kind of two body kinematics such interaction 
does not change fraction of projectile momentum carried by the valence partons. At the same time the distribution 
of small x bare particles is far from understood because the operator of the number of bare particles is not the 
integral of motion in QCD. So, the spectrum of leading hadrons, the probability of the gap in rapidity between 
projectile, and target fragmentation may depend on atomic number. The semiclassical approximation agrees well 
with the conventional picture of high energy hadron-hadron collisions. Fast partons in a nucleon with momentum pjy 
are contracted to a pancake of transverse radius r/v and a small longitudinal size z = r/vmjv /pn while the small x 
partons form a pancake of xy/x times larger longitudinal size (xy ra 0.2 is the average x of the valence quarks). The 
strong interaction of a nucleon with a target occurs when the longitudinal size of the fast nucleon becomes comparable 
with nucleon radius. It is small x partons of longitudinal size ~ 1/m which eventually interact with the target. Thus 
valence quarks and gluons of the projectile keep practically the same longitudinal momentum during collision. It is 
the distribution of small x partons which is influenced by the nucleus medium. This phenomenon is the essence of 
the parton model, and of the Dokshitzer-Gribov-Lipatov-Altarclli-Parisi (DGLAP) evolution of hard processes and 
current models for the energy losses [11]. In the next discussion, we will use this property of QCD. 

A parton with sufficiently large x belonging to a fast proton emits a virtual photon (hard gluon) long before the 
target and it interacts with the target, in a black regime releasing the fluctuation, e.g. a Drell-Yan pair. This leads 
to a qualitative change in the interactions for the partons with x p ,pt satisfying the condition that 

x A = Ap 2 t /(x p s NN ), (1) 
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is in the black body kinematics for the resolution scale pt < p b - bL (xa). Here, (xa) is the maximum p t for 
which the black body approximation is applicable. In the kinematics of LHC, Q 2 « 4(p£ ) 2 can be estimated by 
using formulae derived in [12]. At minimal xa p bbL may reach 4 GeV/c. All the partons with such x p will obtain 
Pt(jet) ~ pf b L (xA) leading to the multijet production. The blackbody regime will extend down in x p with increasing 
the incident energy. For LHC, for p t < 3GeV/c, this regime may cover the whole region of x p > 0.01 where of the 
order of ten partons reside. Hence, in this limit most of the final states will correspond to multiparton collisions. For 
pt < 2GeV/c the region extends to x p > 0.001. At LHC fragmentation of such partons result in the generation of 
hadrons at central rapidities. The dynamics of conversion of the high p t partons with similar rapidities to hadrons is 
certainly a collective effect which deserves a special consideration. 

The total inclusive cross sections in pA scattering can be calculated within BBL using a similar method to that 
used in 7* - nucleus scattering [13]. In particular, the total cross section of the dimuon production is: 

d<j(p + A^^+ f i-+X) Ana 2 K(x A ,x p ,M 2 ) 2 1 2 2 

dx^dx~ p = 1 NP F2p{Xp > Q } ' 6^ M ' 2nRA H*o/*a). ( 2 ) 

Here the if-factor has the same meaning as in the leading twist case, but K — 1 should be smaller since it originates 
from the gluon emissions from the parton belonging to the proton only. x$ is the maximal x for which BBL is valid. 
For the smallest xa (forward kinematics) Eq.2 may be valid at LHC for M 2 < 60GeV 2 . Obviously, Eq.2 prediction 
for the M 2 , xa dependence of the cross section are distinctively different from the DGLAP limit. This difference 
would be less pronounced in the case of pp scattering where scattering at large impact parameters may mask the 
BBL contribution. Another signal for the onset of the BBL is a broadening of the pt distribution of the dimuons as 
compared to the DGLAP expectations, see [14] for a calculation of this effect in the color glass condensate model. 
This effect is similar to the case of p t distribution of leading partons in the deep inelastic scattering [13]. 

As xa decreases further, the formulae for BBL will probably overestimate cross section because the interaction with 
heavy nucleus of sea quarks and gluons in the projectile proton would become black as well. 

The onset of the BBL will lead also to gross changes in the hadron production: there is a much stronger drop with 
xp of the spectrum in the proton fragmentation region accompanied by a significant p t broadening of the spectrum. 
There is also the enhancement of hadron production, at smaller rapidities. Indeed, the individual partons in this limit 
are resolved at the virtuality scale corresponding to transverse momenta <~ p b - bl - without losing a finite fraction of the 
light cone momentum. Hence, the limit of independent parton fragmentation [15] will be realized with an important 
amplification since the leading partons will have much larger transverse momenta [16] than that expected from the 
the estimate of [15] based on the p t broadening observed at the fixed target energies. We want to emphasize here that 
the the approximation of zero fractional losses holds in the pQCD regime and appears to hold in the various models of 
the onset of the BBL, see e.g. [16]. At RHIC this regime may hold for the very forward hadrons and could be checked 
[17] by studying the production of leading hadrons in the central p( 2 H)A collisions [18]. The propagation of a proton 
interacting in the BBL results in the removal of all partons in the nucleus with the proton impact parameter and 
Pt < p\ h1 ' ■ This leads to the formation of a ~ 1 fm radius tube in a perturbative phase. Thus, the essence of the BBL 
is the striping of nucleons off the soft QCD modes and releasing the gas of free quarks and gluons with large transverse 
momenta. Remember that the pQCD interactions between quarks and gluons within the same fragmentation region 
is small. Thus, we conclude that the effects of the BBL should be first manifested in the projectile fragmentation 
region and should gradually expand towards central rapidities. The detectors with a forward acceptance would be 
optimal for this physics. 

Let us now discuss the nucleus fragmentation region in the central heavy ion collisions. It was discussed previously in 
the framework of the soft dynamics, see [19] and references therein, and a significant increase of the densities was found. 
Let us demonstrate now that even more striking effects are expected in the BBL regime. In this case (in difference 
from pA collisions) soft modes will be stripped off in the whole volume of the nucleus. To calculate the properties of 
the quark - gluon state produced in the tube of radius Ra in the target fragmentation region we first evaluate the 
main characteristics of the process in the rest frame of the fragmenting ion. The incoming nucleus is a pancake shape 
with the longitudinal length <~ 1/m for the soft interactions and the longitudinal length r at (to at /pn){xv I 'x) for hard 
interactions. Hence, the nucleons at different locations along the collision axis are hit one after another. In the BBL 
no spectators are left. The hit partons are produced with the same x that they had in the nucleus (since the energy 
losses are oc 1/s). We also have average, p t ~ p bbL and virtuality, p? < (p b - b - L ) 2 . The partons move in the direction 
of the projectile nucleus. Since they are emitted at finite angles their longitudinal velocity is smaller than the speed 
of light. Hence, they are left behind the projectile wave. However, since the emission angles are small the shock wave 
is formed compressing the produced system in the nucleus rest frame, see discussion below. 

To estimate the produced densities, we first calculate the emission angles. Since the parton's x is not changed, 

{Ei- Pi) = xm N , (3) 
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leading to 



p z = (p 2 + p 2 )/2xm N - xm N /2 w p 2 /2xm N . (4) 

Here in the last equation we have neglected /j 2 compared to p\ which is legitimate in the leading order. Since fi 2 > 0, 
neglected terms would increase p z making the emission angles, 9, even smaller. Thus in the BBL the angles 9 

^ Pt/Pz ~ 2xm N /p t , (5) 

are small. So the length of the produced wave package is reduced from a naive value of 2 Ra by a large factor 

S = 1/(1 - cos 9) w p 2 /2x 2 m 2 N . (6) 

However, we must also take into account that the products of the nucleon fragmentation as a whole move forward 
in the target rest frame. The four-momentum of the system can be calculated from the condition that this parton 
system carries almost all of the light-cone momentum of the initial nucleon: 

(^M 2 +p 2 z - Pz )/m N = l, (7) 

where M 2 — J2iPi t/ x i 1S ^ e square of the invariant mass of the system. Hence, p z = M 2 /2rriN, and the Lorentz 
factor 

7 = E/M = y/M 2 + (M 2 /2m N ) 2 /M w M/2m N . (8) 

Combining Eqs.6, 8 we find for the decrease of volume relative to the nucleus volume: 

D = (2m N /M)- < p 2 /2m 2 N x 2 > . (9) 

First, to make simple numerical estimates we assume that all relevant partons carry equal light cone fractions 1/N. 
(The dispersion in Xi leads for a further decrease of the volume and hence to a larger density of the produced system, 
see below.) In this case we obtain 

D = M/m N = Np t /m N . (10) 

Thus, we have demonstrated that the volume is indeed significantly smaller than the nuclear volume. 

The next step is to estimate the energy per unit volume, R E . It is convenient to present it in units the energy 
density of nuclear matter which equals 0.16 GeV/fm 3 . The lower limit on R E is obtained by neglecting the dispersion 
in Xi and taking Xi = 1/N, 

R E = D- N Pt /m N = N 2 p 2 /m 2 N . (11) 

At LHC for x > 0.01 the BBL extends up to p t = 2 3 GeV/c, and the number of partons for such x and virtualities 
is ~ 8 -7- 10. This leads to 



R E > 250. (12) 
If we take into account the dispersion in Xi we find: 

RE = J2Pit/m 2 N xl (13) 

i 

Though in the first approximation the dispersion can be neglected for valence quarks and probably for the gluons 
separately, the quarks carry an average x that is a factor of two larger than gluons. This leads to a further increase 
in the energy density. Taking N q — 3,x q — 1/6, N g — 6,x g = 1/12, we find a rather modest increase of Re (a factor 
of 4/3) for the same N = 9, corresponding to the energy densities > 50GeV/ fm 3 . (If we assume proximity of BBL at 
RHIC for the fragmentation region for p t ~ lGeV/c, wc find quark-gluon energy densities ~ WGeV/ fm 3 .) [20] More 
significantly, the difference in average x's of quarks and gluons leads to a different direction of the flow of the quarks 
and gluons in the cm. frame of the produced system. For the above numerical example, k^/k t ~ 0.6 for quarks and 
<~ —0.3 for gluons. Obviously, this pattern will enhance the interactions of quarks and gluons at the next stage of the 
interactions, making equilibration more likely. 

Let us briefly discuss these interactions and the possible experimental signals. It follows from Eqs.9 - 13 that at 
LHC in the first stage of collisions a strongly compressed hot quark-gluon state of the ellipsoidal shape is formed 
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with the small principal axis of ~ 1/m and density p > 25 partons per /m 3 . At the higher rapidity end, this 
ellipsoid boarders essentially parton free space; on the end close to central rapidities, it boarders a hot qqg state. The 
scattering length for parton i can be estimated as k — 1/ corresponding to the scattering length being 
smaller than 1 fm for <r > 0.5mb. To estimate the interaction cross section, we note that the average invariant energy 
s w 2p1 ~ 8GeV 2 . The initial stage of reinteractions certainly is a highly non-equilibrium process. Nevertheless to do 
a pcrturbative estimate we can conservatively introduce a cutoff on the momentum transfer p ~ ^p -1 / 3 , leading to 
the leading order estimate for the gluon - gluon cross section > lmb. Nonperturbative effects, which remain strong in 
the gluon sector up to y/s ~ 3 GeV, are likely to increase these interactions further. Consequently, we expect partons 
to rescatter strongly at the second stage, though much more detailed modeling is required to find out whether the 
system may reach thermal equilibrium. The large angle rescatterings of partons will lead to production of partons at 
higher rapidities and re-population of the cool region. In particular, two gluons have the right energies to produce 
near threshold cc pairs and in particular x c -mesons with rather small transverse momenta and xf(cc) ~ 2x g ~ 0.1. 
Also leading photons can be produced in the qg — > ^yq subprocesses, though in difference from the central region the 
qq — > p + n~ production will be suppressed due to the lack of antiquarks. Another high density effect is production of 
leading nucleons via recombination of quarks with subsequent escape to the cool region. Hence we expect a rather 
paradoxical situation that the production of leading hadrons in AA collisions will be stronger than in the central pA 
collisions. 

To summarize, we have demonstrated that the onset of the black body regime in the interactions in the target 
fragmentation region which is likely at LHC for a large range of virtualities, will lead to the formation of a new 
supcrdense initial state in the nuclei fragmentation region with densities exceeding nuclear densities at least by a 
factor of 300. Our reasoning is however insufficient to demonstrate whether thermalization processes will be strong 
enough for the system to reach equilibrium necessary for formation of metastable states - new QCD phases suggested 
in a number of recent papers, see the review in [21]. 

Further studies of the experimental manifestations of the formation of high density quark-gluon states are needed 
as well as an investigation of hard processes in the central proton-nucleus collisions which will allow a determination of 
whether this new state of valence quark-gluon matter can be formed at RHIC energies. Similar effects, like correlation 
of high p t hadron production in the fragmentation regions, should be present in the central pp collisions at LHC 
(where jet production at central rapidities is used as a trigger of the centrality [9]). 
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